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°° 1 Abstract 



£> ■ We calculate the penetration depth A in the a, b and c directions for a simple 

in 

model of YBa2Cu307. In this model there are two layers — representing a 



Q ■ Cu02 plane and a CuO chain — per unit cell. There is a BCS-like pairing 

in 

0^ ■ (both s wave and d wave are considered) interaction localised in the CuC>2 



planes. The CuO chains become superconducting at temperatures lower than 



C3 

e 

*^ ■ T c because of their proximity to the planes, and there is an induced gap in 

a: 

the chains. Since the temperature dependence of the penetration depth in 



the b direction (along the chains) is sensitive to the size of the induced gap, 
the difference between the shapes of the penetration depth curves in the a 
and b directions reveals a great deal about the nature of the condensate in 
the chains. We find that in our proximity model there are always regions of 
the chain Fermi surface on which the induced gap is much smaller than T c , 
so that the temperature dependence of A;, is always different than that of A a . 
Experimental observations of the of the ab anisotropy show nearly identical 
temperature dependences. The main result of our paper, then, is that a simple 
proximity model in which the pairing interaction is localized to the planes, 
and the planes are coherently coupled to the chains cannot account for the 
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super fluid on the chains. 
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I. INTRODUCTION 



It is widely believed that the source of the pairing interaction which is responsible for the 
superconducting trasition in the high T c cuprates lies in the CuC>2 planes, which are common 
to all of the cuprates. Many models which attempt to explain high T c superconductivity are 
two-dimensional, which is a reflection of the assumption that the only active pieces of the 
crystal are the CuC>2 planes and that the remaining ions act as placeholders or as charge 
reservoirs. In some materials, however, there are additional layers whose behaviour is not 
clear. In E^S^CaC^Og, for example, it has been suggested that the BiO layer plays the 
role of a normal metal in close proximity to a superconducting material.0 

The only materials in which there is clear evidence that the CuC>2 planes are not the 
only active portion of the unit cell are YBa2Cu3C>7 (Y-123), YBa2Cu40g (Y-124) and their 
close relatives. In these materials there are quasi-one dimensional CuO chain structures. 
Experiments measuring the d.c. resistivity,!'! the infrared and optical conductivity,@ and 
the penetration depth in untwinned crystals!! and ceramics! have found large anisotropies 
between the a direction (in-plane, perpendicular to the chains) and the b direction (in-plane, 
parallel to the chains) which suggest that substantial currents are carried along the chains 
in both the normal and superconducting state. 

In the superconducting state, the source of the condensate on the chains in unclear. 
One possibility is that the pairing interaction is localized to the Cu02 planes, but that the 
chains become superconducting by a proximity effect. In the proximity effect, an intrinsically 
normal metal which is in close contact with a superconductor becomes superconducting near 
the junction as a result of pair tunneling through the junction. The size of the induced gap 
in the normal metal depends on the strength of the coupling across the junction. Y-123 and 
Y-124 are good candidates for proximity effect models because they have the least anisotropy 
between the in-plane and c axis transport properties of the cuprate superconductors, and 
should therefore have a relatively large coupling between the chains and planes. 

Proximity effect models have been studied in the context of high T c materials for a 
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number of years. The most common point of view is that the unusual properities of the 
cuprates can be explained by an isotropic BCS pairing interaction which is contained in 
one of the planes.SEHll The idea behind most of the work is that although the pairing 
interaction may be inherently isotropic, the strongly anisotropic band structure leads to a 
gap structure which may account for the unusual superconducting properties of the cuprates. 
The current authors have taken a different point of view in recent work on proximity effect 
models.EUll We have assumed that the pairing interaction in the planes is intrinsically d 
wave and then attempted to assess the influence of coupling to the chains. Closely related 
to the proximity models are the S/S' multilayer models in which there are two (or more) 
different superconducting layers in the unit cell. There have been detailed examinations of 
the roles played by interplane and intraplane pairingllHl! and a few quantitative calculations 
of physical properties,i@0 but these models have been less thoroughly explored than the 
proximity model because of their relative complexity. 

In this article we address the issue of whether a proximity model can account for the 
condensate on the chains in the YBaCuO compounds. We do this by calculating the pene- 
tration depth A for a simple s or d wave model in the a, b and c directions. In particular, we 
are interested in comparing the temperature dependences of A a and A& with experiment.il In 
our model, the unit cell consists of a Cu02 plane layer and a CuO chain layer. The Cu02 
planes contain the pairing interaction and are coupled to the CuO chains through coherent 
single electron tunneling so that there is an induced gap in the chains. Calculations of the 
penetration depth in a similar model have been made before,0lll although the intrinsically 
normal layers were planes and not chains, and the emphasis was on the anisotropy between 
X ab and A c . 

In Sec. we introduce our model Hamiltonian and find the single particle Green's func- 
tions which we will need for the penetration depth. In Sec. |TTT| we derive an expression for 
the penetration depth which is suitable for a two band, tight binding model. The calculation 
differs slightly from one we made previously^ In Sec. [TV] we discuss the results of numerical 
calculations of the penetration depth, and in Sec. [V] we broaden the scope to a discussion of 



the nature of the condensate on the chains. 



II. HAMILTONIAN 



The goal of this section is to introduce our model for YBCO, and to find the single 
particle Green's functions necessary for the calculation of the penetration depth in Sec. 
IFJ. We begin with a Hamiltonian which describes a system with two layers per unit cell. 
Adjacent layers are separated by a distance d/2. The first layer represents a CuC>2 plane and 
it contains a BCS-like pairing interaction. The second layer represents a CuO chain. It has 
a one dimensional dispersion and is intrinsically normal. The chains are superconducting, 
however, because of their coupling to the planes through single electron tunneling. The 
Hamiltonian, expressed in the Nambu formalism, is: 



H - N/i = C t (k)Q(k)C(k) + const., 



(1) 



where 



Cfk) 



'1-kJ. 



-2kT 



(2) 



Q 



(3) 



and 
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r(k) o 6(k) o 
o -t(-k) o -£ 2 (-k) 

This Hamiltonian has been discussed at length elsewhere0'0El and we only describe it briefly 
here. The dispersions £i and £2 are for the plane and chain layers respectively. We assume 
tight binding dispersions so that £1 = — 2cri[cos(fc a; a)+cos(A; 2/ 6)]— /xi and £2 = — 2o"2 cos(k y b) — 
H21 where a and b are the lattice constants in the planes. In optimally doped Y-123 crystals, 



a and b differ by ~ 1%. For the numerical calculations done in this work, we take a = b. We 
take o"i = 100 meV so that the full bandwidth of the CuC>2 plane is 0.8 eV. For the chain 
layer we take 02 = 80 meV so that \ 2 a {T = 0)/Aj(T = 0) ~ 2.5, as seen experimentally 
by Basov et aM The chemical potentials are /xi = —80 meV and fi2 = 40 meV, which 
yields a Fermi surface that is in qualitative agreement with band structure calculationsJ^H 
We have absorbed an arbitrary band offset into the chemical potentials so that /ii 7^ ^2- 
The chains and planes are coupled by the matrix element t(k z ) = — to cos(/c 2 <i/2), where 
d/2 is the distance between the chains and planes. The chain-plane coupling affects the 
penetration depth in two ways. First, t(k z ) determines the c-axis transport properties. In a 
previous paper0 we have shown that the ratio A C (T = 0)/ X a b(T = 0) varies inversely with to- 
Second, t(k z ) determines the size of the superconducting gap induced on the chains, which 
is reflected in the temperature dependence of Aj,. In this article we choose to = 50 meV 
which yields \ 2 C {T = 0)/X 2 ,(T = 0) ~ 100, which is in rough agreement with experimental 
observations. IHEI 

The final feature of our Hamiltonian is that there is a pairing interaction in the plane 
which drives the superconducting transition. As we have mentioned above, the chains also 
become superconducting at T c through their coupling to the planes. The pairing interaction 
in the planes has the form Vkk' = VVkVk' with rj^ = 1 for an s-wave superconductor and 
r/k = cos(k x ) —cos(k y ) for a d x 2_ y 2 superconductor. Since the pairing interaction is separable, 
the order parameter 



(where Q is the volume of the crystal) can be written Ak = A r/k- 

Diagonalization of the Hamiltonian leads to four energy bands Ei = E + , E 2 = E_, 
E 3 = -E-, E A = -E+ with 




X]^kk'(ci-k'iCik'T) 



k' 



(4) 



Ei 



+ a + A k 



± 



r^-ei+Akf 




(5) 
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In Fig. [l| we show the Fermi surface for a range of k z values between and n/d. The Fermi 
surface consists of two surfaces on which E_ vanishes in the normal state. The two surfaces 
are given by the two solutions to £i(k)£ 2 (k) = t{k z ) 2 . When k z = ir/d, t{k z ) = and the 
two pieces of Fermi surface are those of the isolated chain and plane subsystems. When 
t{k z ) 7^ 0, the chain and plane states form hybrid bands whose energies are given by 



6 + 6 



± 



\ 
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in the normal state. The shift in the band energy due the chain-plane coupling is clearly 
dependent on the relative sizes of t 2 and (6 — 6) 2 - The effect of the chain-plane coupling 
on the Fermi surfaces shown in Fig. [I] is largest in the neighbourhood of the Fermi surface 
crossing at £i = £2 = 0. 

The quasiparticle operators in the diagonalized representation are 

i 



e i (k) = x;t^(k)c i (k), 

where U(k) is the 4x4 matrix which diagonalises Q: U = \U\ U2 U% U4] with 



tA 
tB 



(6) 



(7) 



A = t 2 -(A k + ^ + ei)(^ + 6) 



B = t 2 -(Ai + E J -^)(E 3 -^ 2 ) 



C = A 2 [t 2 + (Ej - 6)1 + B 2 [t 2 + (Ej + 6)1- 



Now that we have diagonalized the Hamiltonian, we can find the single particle Green's 
functions which we will require in the followin section. Defining the temperature Green's 



function G(k; r)^ = — (l/7i)(TCj(k; — rr)Cj(k; 0)), where T is the fermion time-ordered 
product, we have 

m=l l Kl - ^m(KJ 

where 



JO 



(i = (21 + l)n/h(3 are the fermion Matsubara frequencies and (3 = l/k B T, where k B is the 
Boltzmann constant and T is temperature. In our calculation of the penetration depth we 
will need 



G(k;r = 0-),, = i^^(k)/[^(k)]C/i(k), (9) 
Ai i=i 

where f(x) = 1/[1 + exp(fix)]. 

We finish this section with a brief discussion of our usage of the word "gap" . In the model 
presented above, there is only one order parameter, A k , and it describes the condensate in 
the CuC>2 planes. For a separable potential, Ak has the symmetry of the pairing interaction. 
In a multiband material, however, A k is not simply related to the pair wavefunction. For 
example, the anomalous Green's function (which is essentially the pair wavefunction) in the 
Cu02 plane is 

Notice that the symmetry of G± 2 is not the same as the symmetry of A k . For this reason, the 
term "gap" is kept distinct from the term "order parameter", which refers to A k . Perhaps 
the most useful working definition of "gap" is that it is the value of E_ on the Fermi surface. 
Clearly, by this definition, the gap is k-dependent. In regions of the Brillouin zone where 
a section of Fermi surface has predominantly chain (or plane) character, the gap can be 
associated with the chains (or planes). It is wrong to think of the pairs being localised to 
the chains or planes, however; the pairing amplitude between an electron in the chains and 
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an electron in the planes (measured by Gu and G23) is nonzero. In fact the picture of a 
gap belonging to a plane or a chain breaks down in regions of the Brillouin zone where the 
Fermi surface is a strong hybridization of the chain and plane bands. 



III. PENETRATION DEPTH 

The penetration depth is found using an approach which is suitable for the tight binding 
limit. This approach is slightly different than that of our previous work,0 although it yields 
quantitative results which are nearly identical. As we shall see, however, the current method, 
which is based on one used by Hirsch and Marsiglio for a one band tight binding model,ill 
is more satisfying from a physical point of view. 

We begin by writing out the current operator j in the absence of a magnetic field:il 



jo(q = 0) = J d 3 r j( 



f[HS,P], (10) 



where 



P-^RiC^ROc^Ri), (11) 

i,Ri 

is the polarization vector and H n is the Hamiltonian in the normal state. The operator 
cf(Ri) creates an electron in the Wannier state located at the sublattice point Ri. The set 
of points {Ri} refers to the plane sublattice while {R2} refers to the chain sublattice. The 
Wannier representation is connected to the k-space representation by 

c^H^E^c^k), (12) 



N k 



where N is the number of lattice sites. 

In the normal state the Hamiltonian, Eq. (fj), can be written in the Wannier representa- 
tion as: 
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H o = "E E ^cKRi + rOc^RO 

i=l Ri,r; 



- | E H(Ri)c 2 (R 2 

Ri,R 2 
+ "Ri-zd/2,R 2 



R 



+ H.c.} 



(13) 



The vector ri is the displacement to the nearest neighbours of Ri within the plane, z is the 
unit vector in the z-direction, and H.c. indicates the Hermitian conjugate. This Hamilto- 
nian describes nearest neighbour hopping both within and between the chains and planes. 
Substituting Eqs. (0) and ([13]) into Eq. (|10D we get 

jo = J E °i r i c l( R i + ri)ci(Ri) 



i,Ri,ri 



iet d 
~4h 



^zfc^ROc^Ri + zd^) 

Ri 



c t 1 (R 1 )c 2 (R 1 -M/2) - H.c. 



(14) 



In the presence of a finite magnetic vector potential A(r), the tight binding Wannier 
states are modified by a phase so that 



Cj(Ri) — > Ci(Ri) exp 



le 

he 



Ri A(Ri) 



(15) 



where c is the speed of light. The assumption is made that the vector potential is slowly 
varying over the length scale of the crystal lattice, and we will make use of the fact that 
A(q) is strongly peaked about q = throughout this section. To linear order in A, then, 
Eq. ( |T3"D becomes 



H n = Hq — —jo • A(q = 0) 



(16) 



and Eq. (0) becomes 



j=jo-4- E ^ r i[ r i" A(Ri)]ct(R i + r i )c i (R i ) 

he 



i,Ri,ri 



e 2 t d 2 
8h 2 c 



Ez[z- A(Ri)] [c t 1 (R 1 )c 2 (R 1 + zd/2) 



Ri 



+ c t 1 (R 1 )c 2 (R 1 - zd/2) + H.c. 



(17) 
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In the presence of a magnetic field, the observable current is given by (j), and not (jo). We 
can rewrite Eq. (|TTD in a k-space representation using Eq. (O): 



J(q = o) = fE[c t (k)7kC(k) 



ft k 

e -Ct(k)[^-A(q = 0)]C(k) 



(18) 



hcVt 

The vector 7k is a 4 x 4 matrix with three spatial components. It is essentially the Fermi 
velocity: 

where Qij is the Hamiltonian matrix of Eq. (||) and the factor (— comes from the fact 
that Cj annihilates electron states for i = 1, 3 and hole states for z = 2, 4. The dyadic 7 k is 
also a 4 x 4 matrix and it is essentially the effective mass tensor: 



7kij 



dkdk' 

The first term in Eq. ( |i~8D is jo(q = 0) while the second term contains the remaining two 
terms in Eq. (|I7D . 

The current (j) which is generated by the applied magnetic field is given, to linear order 
in A by the Kubo formula:^ 

(i(*)> = (K*)>o + U t dt' ([--jo(t') • A(0, jo(t)\)o. 

n j-oo c 

In the London limit, for the case of a static applied field, this gives 

(j(r)> = 0(r))o-^£GyO,0;0)A y (r), (19) 

where 

1 



hfi Jo 

x (Tjo M (q,-ir)jo iy (q / ,0)) . 

G- 7 is the current -current correlation function, T is the boson time-ordered product, uj\ = 
2lir/h(3 are the boson Matsubara frequencies and // and v refer to the spatial components 
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of jo- The expectation values ( )o are taken with respect to the zero field wavefunction. It 
is straightforward to evaluate Gr 7 in terms of the single particle Green's functions: 



(-> 

G%,(P, 0; iui) = lim — T - ^ Tr [G(k; i( n - iut)j kix 
x G(k + q; i( n )% u ] . 



(20) 



The first term in Eq. (|19[) is the diamagnetic contribution to the screening current. Using 



Eqs. (H) and (|18l) we may evaluate this explicitly: 



1 



j( r )Ua= ^E0(q))^a e 



vth z c 



ETr[G(k;r = 0-)7 k 



AM 



£E/[^(k)]7L-A(r) 



(21a) 



k i=l 



where Tr is a trace over the components of the 4x4 matrix contained in the square brackets, 
and 7 k = t/^(k)7 k £/(k). In order to derive Eq. ( |21a| ), we have used the fact that (jo)o = 
and that A(q) is peaked about q = 0. The second term in Eq. (|i~9"D is the paramagnetic 
contribution to the screening current. Evaluating Eq. (|20"D explicitly we have 

„2 4 



(j) 



para 



-^2~ c E EW7kji- A(r)] 



X 



i,j=l k 

df(Ei) 
BE,. 



[1-5. 



1,31 



J{ E i )-f{E j 
Ei — E~ 



(21b) 



where 7k = £/t(k)7k?7(k). This expression for the paramagnetic current is the same as in 



Ref. [LJ where it was discussed at length. We will only repeat the points which are directly 
relevant to the current work, and the interested reader is referred to our earlier work. The 
total current produced by the magnetic field is 



(j) = (j)dia + (fipara 
= -K A. 



(22) 
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It is straightforward to show that = if // ^ v (recall that /i and v refer to spatial 
directions) so that the penetration depth is given by 

4 = *W (23) 

This is the main result for this section. In order to plot A~ 2 as a function of temperature in 
Sec. [IV], we must evaluate the integrals in Eqs. (|21a| ) and (|21b| ) numerically. 



We will finish this section with a few comments about Eqs. ( |21a| ) and (|2 lb|) . In the usual 
treatment of the penetration depth the diamagnetic contribution to the screening current is 
$ dia = — ne2 A/mc, which is independent of temperature. The temperature dependence of 
the penetration depth, then, comes from the paramagnectic contribution to the screening 
currents which, for a one band free electron metal with an isotropic gap, is 

where E — [e 2 + A 2 ] 1 / 2 , is the Fermi velocity, and fj, is the chemical potential. The 
paramagnetic term counts the number of thermal excitations (broken Cooper pairs) which 
degrade the screening current. At T = the paramagnetic term vanishes, so that j = 
—ne 2 A/mc. When A = the paramagnetic term cancels the diamagnetic term exactly so 
that j = 0. For systems which are more complicated than the free electon gas, it is common 
to make the approximation = —jpara A ■ The approximation is exact at T = T c and, 
provided the temperature dependence of the diamagnetic term is weak, the approximation 
is a good one. This is the approximation we made in our previous discussion of the two 
layer model.B In the current work, however, we have treated Sdia * n a f asn ion which is 
more consistent with the tight binding model, so that while Eq. (|21b|) is the same as we 
found previously, Eq. ( [2 lap is different. There is little quantitative difference between the 
two approaches, however, since both expressions for the diamagnetic current are weakly 
temperature dependent and both cancel the paramagnetic current above T c . 

The most significant difference between Eq. Q21ap and the usual expression for the pen- 
etration depth is the interband term, which is proportional to \f{E,j) — f(Ej)]/[Ei — Ej]. 
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While the intraband term (which is proportional to df(Ei)fdEi) counts the number of ther- 
mally broken pairs, the interband term describes the degradation of the screening currents 
by interband transitions. The interband term does not vanish at T = so that, unlike the 
single band case, there is a finite paramagnetic contribution to the screening current. 



IV. RESULTS 

The question we are attempting to address in this article is whether a proximity effect 
model can account for the experimentally observed anisotropy in the temperature depen- 
dence of the penetration depth in Y-123. In this section we will present the results of 
numerical calculations of the penetration depth for the model Hamiltonian introduced in 
Sec. 0. We will compare these results to experiments and to related calculations made 
with a two-plane proximity model (in which the intrinsically normal layer is a two dimen- 
sional plane). One of the main goals of this section is to emphasize the difference between 
two-plane proximity models and chain-plane models of the type studied here. 

To begin with, we will discuss calculations of the penetration depth in the two-plane 
proximity models.00 One of the important features of the proximity model is that it in- 
troduces low energy excitations into the superconducting spectrum. The reason for this is 
that the induced gap in the intrinsically normal plane is proportional to the strength of 
the chain-plane coupling t(k z ) (which vanishes at k z = n/d) so that the gap will have a 
nodal structure even if the pairing interaction has isotropic s wave symmetry. The need 
for a gap structure with nodes has been suggested, for example, by measurements^! of X a b 
(the in-plane penetration depth) in twinned single crystals of Y-123. The linear dependence 
of X a b{T) on T at low temperatures is easily explained by any gap structure with nodes!! 
While these measurements are commonly taken as support for d wave models,^ it has also 
been shownffl that two-plane proximity models also result in linear low T behaviour. Since 
a central theme in much of the work on proximity mo delJHoffi is that the pairing in- 
teraction in the intrinsically superconducting plane is s wave, the low energy excitations in 
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the induced gap are an essential feature of the proximity models. 

In Fig. 0(a) we plot the penetration depth for our plane-chain proximity model for 
the case of an s wave gap. We find that, unlike the case of the two-plane model, the 
temperature dependences of A a and Xb are dramatically different. The most important 
difference is that the temperature dependence of A a is nearly identical to that of a single 
layer s wave material with no chains, while Xb has a linear low T behaviour similar to that 
found in the two-plane proximity models.0 The factor of two difference between A a (0)~ 2 
and Af>(0)~ 2 comes from the screening currents carried in the b direction by the chains, and 
the linear T dependence in A^ 2 at low temperatures comes from the node in the induced 
gap at k z = n/d. The fact that the low T behaviour of A~ 2 is exponential and not linear 
indicates that pairs associated with a axis screening currents have a finite gap for all values 
of k z . We can understand this in more concrete terms as follows: In Sec. [Ill] we showed 
that the screening current has two parts — a diamagnetic part which is roughly independent 
of T and a paramagnetic part which accounts for processes (such as thermal pair breaking) 
which degrade the screening currents. The temperature dependence of the penetration depth 
comes from the paramagnetic screening current, given in Eq. ( |21b[ ). Despite its complicated 



appearance, Eq. ( |21bp has a simple physical interpretation. The factors 7k are electron 



Fermi velocity vectors, while the two terms involving Fermi functions count the number of 
thermally excited quasiparticles which participate in intraband (i = j) or interband (i ^ j) 
paramagnetic processes. When we calculate the screening current in the a direction, then, 
the integrand in Eq. (|21b| ) is weighted by the square of the Fermi velocity in the a direction. 
In Fig. |l|, we can see that this is small both on segments of the Fermi surface associated with 
the chains and on segments of the plane Fermi surface which are distorted by the chains. 
The most obvious consequence of this is that the chains do not participate significantly in 
carrying currents in the a direction. A more subtle result is that, even though there is a 
node in the induced gap in the chains, it is not seen by electrons travelling in the a direction 
so that Cooper pairs which are part of the a axis screening current have a finite gap. The 
onset of thermal pair breaking, then, occurs at a much lower temperature in the b axis 
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supercurrent than in the a axis supercurrent. 

In Fig. 0(b) we plot the penetration depth for a d wave order parameter and find results 
which are similar to the s wave case: A a (T) is essentially the same as found in single layer 
d wave models and Af>(T) resembles X a b{T) found in the two-plane proximity models. As 
for the case of an s wave gap, the reason is that there are a larger number of low energy 
excitations in the chains than in the planes. The d wave gap in the planes has nodes along 
k x = ±k y , while the induced gap has nodes along k x = ±k y and k z = ir/d. 

The large temperature dependence of the ab anisotropy seen in Figs. @(a) and (b) is 
difficult to reconcile with measurements of A B (T) and Afc(T) in untwinned crystals. In these 
experiments!! A a and A& have a nearly identical temperature dependence, although their 
absolute magnitude differs by a factor of 1.5 at T = 0. In our model, the temperature 
dependence of the anisotropy is a result of the fact that Cooper pairs in the chains are 
more easily broken than Cooper pairs in the planes. Clearly, then, in a realistic model, the 
density of low energy excitations in the chains must be similar to that in the planes. This 
is not a trivial requirement. It implies that both the nodal structure and the magnitude of 
the gaps in the chains and planes be similar. It is possible to eliminate the nodes in the 
induced gap at k z = njd by, for example, making the ansatz that t{k z ) = t (this would 
describe a single bi-layer). However, this is not sufficient to eliminate the temperature 
dependence of the anisotropy. For regions of the chain Fermi surface where ^> |i(fc z )|, 
the induced gap is of the order0 A k t(/c 2 ) 2 /£ 2 . In Fig. |], the smallest induced gap occurs at 
the intersection of the chain Fermi surface with the Brillouin zone boundary (at k x = it /a) 
at which to/£i(k) 2 ~ °- 02 3. The onset of thermal pair breaking in the chains, therefore, will 
occur at a much lower temperature than in the planes. 

The penetration depth in the c direction as a function of temperature is shown in Fig. |3|. 
The shapes of the curves are similar to what we found in previous work0 in which we exam- 
ined a model with two planes per unit cell. Experimental observations of A c in Y-123@'0il 
and Y-1240 are contradictory. All of the experiments find that at low temperatures A c 2 (T) 
can be fitted by a linear T dependence, A c (0) 2 /A c (T) 2 ~ 1 — aT/T c , but the slope of the fit 

16 



varies dramatically. Two of the infrared experiments!!'!! find that a < 1, while the third§ 
finds that a, ~ 1 and the microwave experiment!! finds that a ^> 1. Until some sort of 
consensus is achieved, it will be difficult to say anything about our model. 



V. CONCLUSIONS 

It is clear that our proximity effect model cannot describe the temperature dependence 
of the ah anisotropy of the penetration depth which has been observed experimentally by 
Zhang et a/.i. Essentially, the problem with our proximity model is that, unlike the case of 
the two-plane model, there are always regions of the chain Fermi surface on which the gap is 
small, so that the temperature scale over which the penetration depth parallel to the chains 
varies is much lower than the scale perpendicular to the chains. The question that needs 
to be answered, then, is to what extent is our model representative of proximity models in 
general. 

The common feature of proximity models is that the chains are intrinsically normal but 
driven superconducting by their coupling to the planes. Where proximity models differ is in 
the nature of the chain-plane coupling. In our model we have made the assumption that the 
chain-plane coupling is coherent, so that chain states are coupled to plane states with the 
same value of k. The amount of mixing between the two states depends on the difference 
in energy between them so that, for example, in Fig. [l| the chain and plane Fermi surfaces 
are most strongly mixed in the neighbourhood of their crossing. In a similar fashion, the 
induced gap on the chain is small (of the order of a few percent of the intrinsic gap in the 
plane) wherever the chain and plane Fermi surfaces are far apart. This is the reason for the 
large difference in the temperature dependence of A a and 

One solution to this is to couple the chains and planes incoherently, so that every state k 
on the planes is coupled equally to every state k' on the chains. There is some evidence that 
there is incoherence along the c axis: Kleiner and Mullerll have found an intrinsic Josephson 
effect in E^S^CaC^Og and, more recently, in underdoped Y-123. The d.c. resistivity of 
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Y-123 in the normal statea shows semiconducting behaviour in underdoped samples, and the 
optical conductivity along the c-axis (see, eg., Ref. |27] and references contained therein) has 
a non-metallic response. For an incoherent model of the type described above, the induced 
gap on the chains is proportional to the average over the Fermi surface of the gap on the 
planes. The difficulty with this model is that, for a d wave order parameter, the induced 
gap in the chains will vanish. It is, in fact, a general feature of d wave order parameters 
that they do not contribute to incoherent processes (see, for example, Refs. ^ and |l^). If 
on the other hand, we assume that the order parameter in the planes has an isotropic s 
wave symmetry, then the induced gap on the chains will not vanish. The problem now, 
however, is that incoherent coupling does not introduce a nodal structure into the gap the 
way coherent coupling does so that it is difficult to reconcile such a model with a linear 
low temperature penetration depth. For a model with incoherent chain-plane coupling to 
successfully describe the low temperature penetration depth, it would have to have an order 
parameter with nodes on the Fermi surface but whose Fermi surface average was nonzero, 
and the induced gap in the chains would have to be of the order of T c so that the temperature 
dependence of A a and A& would be similar. 

Leaving, for a moment, the discussion of the nature of the chain-plane coupling, we will 
now turn to a more conceptual problem — that of the size of the chain-plane coupling. The 
coupling strength t is chosen to account both for the fact that A C (T = 0)/A a (T = 0) ~ 
10EEHI and for the size of the induced gap in the chains. As is well known, the chain-plane 
coupling can degrade T c substantially. We find that for t = 50 meV, T c is only 65% of its 
value at t = 0. It is also difficult to reconcile the picture of weakly coupled two dimensional 
planes with such large values of the chain-plane coupling. In our model the ratio of the 
electron hopping strengths along the c and in-plane directions is to/2ai = 0.25, so that it 
is difficult to imagine that the c axis coupling is a weak perturbation in an otherwise two 
dimensional system. The challenge, therefore, for theories which begin with models of a 
single CuC>2 plane is to explain the large anisotropy between the a and b supercurrents in 
YBCO without invoking a large chain-plane coupling. 
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Kresin and Wol£3 have suggested that proximity effect models require an inelastic chan- 
nel for the chain-plane coupling. In their two-plane model, electrons can hop between the 
planes through coherent tunneling or through scattering from a phonon. Their model is 
more three dimensional than the ones discussed above since the inelastic interplane coupling 
acts as a pairing process which leads to an increase in T c . It is possible that in a chain- 
plane model, some kind of inelastic transport mechanism along the c axis might lead to a 
sufficiently large gap in the chains that A a and A& would have similar T dependences. The 
idea of a mixture of pairing interactions has recently been proposed by Song and Annett,i! 
although they have limited their discussion to mixing phonons and Coulomb interactions 
within a single plane. 

There is also the issue of whether a simple two-band model can be representative of 
Y-123. More careful band structure calculations0'@ find that the Fermi surface has four 
pieces instead of two. The two additional pieces of Fermi surface come from the internal 
structure of the Cu02 bilayer (which we have treated as a single layer) and from the internal 
structure of the CuO chains. The inclusion of these two pieces of Fermi surface is not likely 
to affect the important results contained within this paper however: the additional piece of 
Fermi surface due to the Cu02 bilayer has a nearly tetragonal symmetry (and will therefore 
contribute to the anisotropy in the penetration depth) and the piece due to the CuO chains 
is small and will only make a small change to the screening currents. 

Our final conclusion, then, is as follows: A proximity model for Y-123 in which the su- 
perconducting pairing interaction is localised to the planes and the chain-plane coupling is 
coherent will not account for the temperature dependence of the anisotropy of the penetra- 
tion depth seen in experiments. It is possible that other models for the chain-plane coupling 
will be able to adequately describe the ab anisotropy. The single largest problem faced by 
proximity models is that penetration depth experimentsi seem to indicate that the gap in 
the chains is of the same order as the gap in the planes. 
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FIGURES 

FIG. 1. The Fermi surface for the model Hamiltonian is shown for a range of k z between k z = 
and k z = ir/d. When k z = ir/d, the chain-plane coupling vanishes and the two pieces of Fermi 
surface are those of the isolated chains and planes. As the chain-plane coupling increases, the 
Fermi surfaces hybridize and are pushed apart. The effect of the chain-plane coupling is largest 
where the two Fermi surfaces are closest together. There is an induced gap on the chains whose 
size is greatest where there is the most chain-plane mixing. 

FIG. 2. (a) In-plane penetration depth for an s wave order parameter. The penetration depth 
in the a direction (perpendicular to the chains) is nearly that of a pure s wave superconductor in 
the absence of chains. The penetration depth in the b direction has a very different temperature 
dependence from that in the a direction because the size of the induced gap in the chains is 
much different from the size of the gap in the planes. The relative bandwidths of the chains and 
planes were determined by setting A^(0)/A5(0) ~ 2.5, in accordance with experiment, (b) In-plane 
penetration depth for a d wave order parameter. Again, A a (T) is essentially the same as for a 
single layer d wave superconductor, while the shape of Aft(T) reflects the structure of the induced 
gap in the chains as well as the planes. 

FIG. 3. Penetration depth in the c direction for both an s wave (solid line) and a d wave (dashed 
line) order parameter. The strength of the chain-plane coupling is chosen to be to = 50 meV so 
that Ag(0)/A^(0) ~ 100, as observed experimentally. 
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